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We develop a velocity-dependent one-scale model for the evolution of domain wall networks in flat 
expanding or collapsing homogeneous and isotropic universes with an arbitrary number of spatial 
dimensions, finding the corresponding scaling laws in frictionless and friction dominated regimes. We 
also determine the allowed range of values of the curvature parameter and the expansion exponent 
for which a linear scaling solution is possible in the frictionless regime. 
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I. INTRODUCTION 

Domain walls can form at a spontaneous phase tran- 
sition with discrete symmetry breaking in the early uni- 
verse [ij . Their evolution after the phase transition 0, Q 
plays an important part in determining their possible cos- 
mic role. Although, current observational constraints on 
the dark energy equation of state parameter strongly dis- 
favor domain walls as the single dark energy component 
0i such limits, on their own, do not yet rule out a 
substantial impact of a frustrated domain wall network 
on the acceleration of the Universe around the present 
time. 

The dynamics and cosmological consequences of do- 
main wall networks in flat, expanding or collapsing ho- 
mogeneous and isotropic universes with 3-1-1 dimen- 
sions have been extensively studied in the literature, us- 
ing both semi-analytical models and domain wall 
network simulations [lol - [l^ considering complex models 
with Y- and/or X-type junctions [l5l - [l7| . These studies, 
in combination with recent observational results, have 
provided very strong evidence that domain walls cannot 
account for a significant fraction of the dark energy. 

Until the present date, analytical studies of domain 
wall dynamics in more than 3 spatial dimensions have 
only considered the evolution of maximally symmetric 
configurations fisj . Although a number of interesting re- 
sults have been obtained, they cannot be directly applied 
to the evolution of domain wall networks in homogeneous 
and isotropic universes with an arbitrary number of spa- 
tial dimensions. In this paper we eliminate this short- 
coming by generalizing the existing velocity-dependent 
one-scale (VOS) model, describing the evolution of the 
characteristic length and velocity of a domain wall net- 
work, to an arbitrary number of spatial dimensions. 

The paper is organized as follows. In Sec. II we 
develop a VOS model for the evolution of {N — 1)- 
branes in (N -\- l)-dimensional Friedmann- Robertson- 
Walker (FRW) universes. In Sec. Ill we study the fric- 
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tionless regime, finding a number of scaling solutions de- 
scribing the dynamics of domain walls in in (N -\- 1)- 
dimensional expanding or collapsing FRW universes. In 
section IV we investigate the friction dominated regimes. 
We then conclude in Sec. V. Throughout the work, we 
will assume the metric signature [-|-, —] and the cal- 
culations will be done using units in which c = h = 1. 



II. VELOCITY-DEPENDENT ONE SCALE 
MODEL FOR (iV 1)-BRANES 

In a (A^ + l)-dimensional FRW Universe, the dynamics 
of a featureless {N — l)-brane whose thickness is much 
smaller than its curvature radii is described by Q 



(1 - v^) [NHv - k]=0, 



(1) 



where a dot represents a derivative with respect to phys- 
ical time t, a is the scale factor, H = a/ a is the Hubble 
parameter, v is the microscopic velocity of the surface and 
K is the curvature of the surface defined at each point as 
the sum of the principal curvatures. 

Now consider the case of a network of {N — l)-branes 
(domain walls) in a {N + l)-dimensional FRW Universe 
and let us start by defining the r.m.s velocity v — ^ (u^) 
as 



J v^pdV 



(2) 



where p is the domain wall energy density and V is the 
physical volume. Let us also define the characteristic 
length, L, of the domain wall network as 



P ■ 



(3) 



where (Jat-i is the wall energy per unit — 1-dimensional 
area and p ^ V^^ J pdV is the average domain wall den- 
sity. 

Multiplying Eq. ([1]) by u, making the above weighted 
volume average and then dividing by v one obtains 
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v + -{v{l- v^) [NHv - k]) = . 



(4) 
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If we assume that 
fies to 



then Eq. ^ further simph- 
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= 0, 



where i,,^ — NH, k — kL and 



_ / v{l-v'^)KpdV 
u(l-u2) ~ v{l-v'^) J pdV 



(5) 



(6) 



The above assumption is vahd in the relativistic regime 
up to first order in (1 — f ) and it has a negligible impact in 
the non-relativistic limit. Although it is possible to con- 
struct network configurations with the same v but difl^er- 
ent K, in most physically realistic situations it is sufiicient 
to consider that R — R{v) [T§|. In the presence of fric- 
tional forces the characteristic damping length becomes 
NH + l^^, where if is the friction lengthscale. 
Let us assume that the domain wall network is statis- 
cally homogeneous and isotropic on large enough scales, 
so that it behaves as a brane gas. Energy-momentum 
conservation in a FRW universe implies that 



p + NH(p + P) 



0. 



(7) 



where P = V^^ j V dV and p ~ V^^ J pdV are the av- 
erage domain wall pressure and density whose ratio is 
given by the equation of state parameter [2Q| 



—dt) at = 0. We shall label the various models as 
M/. Here, s — ± depending on whether dt = ztdt^, and 
i = 1,2 or 3, depending on whether /3<0, 0</3<l 
or (3 > 1, respectively. We shall also assume that c > 
0. The models M^, and represent expanding 
solutions with = either at the the big-bang (M^^ 
and M;l) or at the big rip (for Mf ). The models Mi+, 
and Mj^ represent collapsing universes with i» = 
either at the the big-crunch (M^^ and M^) or at an 
initial infinite density singularity with a(t* = 0) = oo (for 
M^). In an expanding universe, if > —1, the Hubble 
radius, increases with time but for —1 < w < Wc 

(with Wc — {2 — N)/N, so that /? > 1) the comoving 
Hubble radius {\H\~^ /a) decreases with time. Note that, 
in an expanding universe, the comoving Hubble radius 
will monotously increase or decrease with time depending 
on whether d is negative or positive, respectively. Of 
course, the reverse is true in a collapsing regime. 
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Eq. ([7]) can be generalized in order to account for addi- 
tional energy loss mechanisms, due to interface collapse 
and friction 



NH (p + V) = 



cv 

T 



p- 



(9) 



Here, c is a phenomenological parameter which depends 
on the specific properties of the network and can be cal- 
ibrated using numerical simulations [13]. Using Eq. © 
one obtains 



L 



(10) 



Eqs. (O and (jlOp constitute a generalization of the VOS 
model in Q to the case of a + 1-dimensional FRW 
Universe, which has been derived here in detail. 



III. FRICTIONLESS REGIME 

For simplicity, we shall assume that the dynamics of 
the universe is driven by a fluid with w = constant 7^—1 
so that a (xt^, where /3 = 2/{N(w + 1)) and > is the 
time elapsed since the initial singularity (if dt* = dt) or 
the time remaining up to the final singularity (if dt* = 



TABLE I: Summary of the main properties of the different M' 
models. The label i takes the value i = 1,2 or 3 for /? < 0, 
0</3<lor/3>l, respectively. On the other hand, s = ± 
depending on whether dt — ±dt, and t, is the time elapsed 
since the initial singularity (s = +) or the time remaining 
until the final singularity (s — —). The remaining -I- and — 
indicate the sign of the cosmological parameters represented 
in the table and the letters C and E represent collapse and 
expansion, respectively. 



A. Linear scaling solutions (£f — 00) 

If the friction lengthscale becomes negligible compared 
to the Hubble radius then Eqs. (O and pUj) may have 
a linear attractor solution. This attractor solution corre- 
sponds to a linear scaling regime of the form 



so that 



L = ^t,, 

(l-/?)fc 
■ Nl3{k + £) 



V = const . 
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FIG. 1: The range of values of the curvature parameter, fc, 
for which the VOS equations admit a Unear scahng solution 
(grey area) , as a function of/3/(l — /3), for the models (left 
panel) and M'^ (right panel), with c = 0.5. The dash-dotted 
(blue), dashed (red) and solid (purple) lines are defined by 
L — du; v'^ = 1/2 and — \, respectively. 



FIG. 2: The range of values of the curvature parameter, k, for 
which the VOS equations admit a linear scaling solution, as a 
function of /3/(l — /3), for the model (grey area). The left 
panel represents the allowed range of /N with c = and 
the right panel shows how this range would be changed for 
c = 0.5. The dash-dotted (blue) and dashed (red) are defined 
by L = dff and = u^ax = 1 ~ V-'^) respectively 



Two necessary conditions for a linear scaling solution of 
the VOS equations to be possible are 



< u < 1, C > 0, 



(14) 



but these are by no means sufficient. There are a number 
of complementary constraints which significantly reduce 
the range of parameters consistent with a linear scaling 
solution. 

The r.m.s. velocity of maximally symmetric {N — 1)- 
branes with a SN~i-i ® II* topology oscillating periodi- 
cally in a Minkowski spacetime is given by [l8| 



N 

N ~i 



(15) 



The value of takes the maximum value, equal to 
1 — l/N, if i = (i.e., all the principal curvatures of 
the surface are equal and non-zero) and the minimum 
value, equal to 1/2, ii i — N — 2 (i.e., if only one of the 
principal curvatures is non-zero). Causality constraints 
do not allow for infinite flat branes to be formed in real- 
istic cosmological scenarios and consequently we did not 
consider the case with i = N — 1, where all the principal 
curvatures are equal to zero, in the above discussion. 

We may use this result to estimate the minimum and 
maximum r.m.s. velocity of a domain wall network in a 
(iV-f l)-dimensional Minkowski spacetime (corresponding 
to /3 = 0) as 
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Note that for /3 — the curvature parameter k must be 
equal to zero for a linear scaling solution with v < 1 to 
be attained. The expansion (collapse) of the universe 
will add a damping (forcing) which will necessarily lead 
to a smaller (larger) r.m.s. velocity, v and a curvature 
parameter, k, larger (smaller) than zero. 



Also, causality constraints require the characteristic 
length of the network to be smaller than the particle hori- 
zon, dn, at any given time 



L < dn = 



dt' 



(17) 



where ti — Q ov ti = —oo depending on whether s = -I- 
or — , respectively. The particle horizon is infinite in the 



case of the models , 



and . The constraint 



given by Eq. (fTT)) is only relevant in the case of the mod- 
els for which a linear scaling solution is possible {M^ , 
M2 and M^) and can be written as ^ < |1 — or, 
alternatively, as v'^ < [k + c)~^. 

The r.m.s. velocity constraints for the models Afj^, 



and M3 determine the range of values of the curva- 



ture parameter, /c, for which linear scaling solutions are 
allowed, for any given /? and N . Figs. [T] and [2] ilhistrate 
this in the case of contracting {M^ and A/3" ) and expand- 
ing {M^) models, respectively. The first thing to notice 
is that, in contracting models, linear scaling regimes are 
strictly forbidden if c = {v^ < is not allowed) . On the 
other hand, if c 7^ the network may attain linear scaling 
regimes with velocities within the physically significant 
range, \ < v"^ < 1. The left and right panels of Fig. [T] 
show the allowed range of the curvature parameter, as a 
function of the expansion exponent, for which the linear 
scaling solutions are permitted, in the case of the models 
(16) and Afj^, respectively (for c = 0.5). In both cases, 

scaling solutions are allowed for every value of /? but the 
allowed range of k is strongly restricted. Also, in both 
models for 2 = 0.5 the causality constraint, given by Eq. 
P?)) . does not introduce further restrictions on the plane 
(/c^/7V, /3/(l — /3)). However, for larger values of c, and 
consequently of ^, the region for which causality is vio- 
lated widdens. In fact, if the value c is sufficiently large. 
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all linear scaling solutions may be forbidden. In the case 
of model , the network is able to reach a linear scal- 
ing regime for j3 > N^^ (or equivalently v < fmax), even 
if c = 0, as the left panel of Fig. [5] illustrates. For larger 
values of /?, the allowed range of the curvature parameter 
is only limited by causality. On the other hand, if c ^ 
the linear scaling regime can be attained for any value of 
/3, as shown in the right panel of Fig. 



B. Inflation and superinflation 



In the case of accelerated expansion, as in the models 
(which is inflationary, with a > 0) and Mj~ (which 
exhibits superinflation with H > 0), a stretching regime 
could persist (note that these models would require that 
c < in order for a linear scaling solution to be possible) . 
However, in that case, the expansion is fast enough to 
hamper the brane velocities and make them arbitrarily 
small. For both models one has 



L cx a , 



0. 



(18) 



In a {N + l)-dimensional FRW universe, in which the 
branes are the dominant component of the energy den- 
sity, the Einstein equations imply 



SttG 



N+l 



N{N - 1] 



N{N - 1) 



{{N~2)p + Np), (19) 
(20) 



where Gn+i is the iV + 1 dimensional Newton constant. 
From Eqs. ([7][H1) one obtains 



N{l + Wb) 1 + 



(21) 



In order to accelerate the universe one needs /3 > 1 or 
equivalently wt < Wc ~ {2 — N) /N and consequently, in 
a domain wall dominated universe, the r.m.s. velocity of 
the walls, 



-2 1 



(22) 



has to be small enough for the universe to be accelerating. 



C. Ultra-relativistic collapsing solution 

Consider the case of model M2 which represents a col- 
lapsing universe with fc < and < /3 < 1. A linear scal- 
ing solution would only be possible if fc < — c. However, 
such solution would necessarily be transient because, if 
locally the curvature scale of the domain walls becomes 
smaller than then they will tend to become frozen 

in comoving coordinates while traveling at the speed of 
light (note that the comoving Hubble radius, \H\^^, de- 
creases with time in this model). 



Though, in general, in the context of VOS models, the 
correlation length, L, may be identified with the typical 
physical distance travelled by a brane segment before en- 
countering another segment of the same size, this identi- 
fication breaks down in the ultrarelativistic limit ^21, . .22i] . 
In this limit, instead, i is a measure of the energy of the 
branes and, therefore, the Lorentz factor 7 = {l — v'^)~^^'^ 
should be included in the definition of the physical length, 
Lph ^ jL. The fraction of the energy lost by the net- 
work due to interface collapse in a timescale dt may be 
estimated as 



dp 
P 



dL 
T 



Lph 



-dt 



-dt. 



Hence, in the u — > 1 limit 

c cx cx . 

As a consequence, we find that 

L^{1 + N)HL, 



(23) 



(24) 



(25) 



which indicates that the brane network will behave effec- 
tivelly as a radiation component. In this limit the domain 
walls are ultra-relativistic with 7 oc but the charac- 
teristic length of the network becomes frozen in comoving 
coordinates so that Lph cx a. The combination of these 
two effects gives the solution to Eq. L cx a^+^. 



IV. FRICTION DOMINATED REGIME 

If a domain wall moves through a radiation fluid, the 
interaction with the ultrarelativistic particles results in a 
frictional force which may be written as [H 



7V 



(26) 



where if ~ o j prad oc a^^^ is the friction lengthscale and 
Prad oc 0"'^+^' is the energy density of the relativistic 
particles. 



A. Expanding universe 

As in the case of cosmic strings and domain wall net- 
works in a (3-|-l)-dimensional universe, we expect (A'^— 1)- 
brane networks in a (A^ -f- l)-dimensional universe to ad- 
mit two different scaling regimes during the initial part 
of its evolution, while the friction lengthscale is signifi- 
cantly smaller than the Hubble radius (£f <C H^^). In 
this limit \v\ ^ v/if so that v = k£[/L. 

If the initial density is sufficiently low {HL ^ v'^L/li 
and HL ^ cv) the branes will be confornially stretched 
by expansion. During this stretching regime the scaling 
laws are 



r «f N 

L CX a , V (X — oca . 

a 



(27) 
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This is a transient regime since the velocity increases 
rather quickly, due to the effect of the domain wall cur- 
vature. 

As the friction lengthscale approaches the character- 
istic length of the network, the Kibble regime emerges. 
Given the fact that, during this regime, the characteristic 
velocity is higher than it has hitherto been, a considerable 
amount of energy will be lost due to brane self-interaction 
{HL ~ cv). Therefore, even though friction still domi- 
nates the dynamics, the scaling laws are different from 
those of the stretching regime 



LaJ^, v^.MH\- (28) 

If the initial density of the network is high enough, the 
Kibble regime occurs right after the formation of the 
domain walls and the network will not experience the 
stretching regime. As the universe expands, the friction 
lengthscale grows faster than (if -I- 1 > and 
will, eventually, overcome the characteristic length of the 
network, thus bringing the Kibble regime to an end. 

B. Collapsing universe 

The scaling solutions given by Eqs. (I?fl) and (^5]) also 
account for the dynamics of a domain wall network in 
a collapsing universe during a friction dominated era. 



In this case, the domain wall network ends in a friction 
domination era, coming to a standstill in comoving co- 
ordinates and then being conformally contracted (with 
L oa a). Hence, in this regime the average energy density 
of the network is given by p oc a~^. As the background 
temperature and density approach those of the original 
wall-forming phase transition, the branes effectively dis- 
solve back into the high density radiation background. 



V. CONCLUSIONS 

In this paper we have studied the evolution of domain 
wall networks, in flat expanding or collapsing homoge- 
neous and isotropic backgrounds with an arbitrary num- 
ber of spatial dimensions, using a VOS model. We have 
obtained the scaling laws in frictionless and friction domi- 
nated regimes and determined the constraints which have 
to be satisfied by linear scaling domain wall networks. 

The present work is a significant improvement over pre- 
vious analytical studies of domain wall network evolu- 
tion, unifying in a common framework the dynamics of 
domain wall networks in expanding/collapsing and fric- 
tionless/friction dominated regimes. The generalization 
of the VOS model to (A^ -I- l)-dimensional FRW universes 
also provides an important tool to describe the evolution 
of domain wall networks in more than 3 spatial dimen- 
sions which, up to now was restricted to very special 
configurations. 



[1] A. Vilenkin and E. P. S. Shellard, Cosmic Strings and 

Other Topological Defects (Cambridge University Press, 

Cambridge, 1994). 
[2] W. H. Press, B. S. Ryden, and D. N. Spergel, Astrophys. 

J. 347, 590 (1989). 
[3] L. Sousa and P. P. Avelino, Phys. Rev. D81, 087305 

(2010). 

[4] J. Frieman, M. Turner, and D. Huterer, Ann. Rev. As- 

tron. Astrophys. 46, 385 (2008), 0803.0982. 
[5] E. Komatsu et al. (2010), 1001.4538. 
[6] M. Bucher and D. N. Spergel, Phys. Rev. D60, 043505 

(1999), astro-ph/9812022. 
[7] P. P. Avelino, C. J. A. P. Martins, and J. C. 

R. E. Oliveira, Phys. Rev. D72, 083506 (2005), hep- 

ph/0507272. 

[8] P. P. Avelino, C. J. A. P. Martins, J. Menezes, 
R. Menezes, and J. C. R. E. Oliveira, Phys. Rev. D73, 

123519 (2006), astro-ph/0602540. 

[9] L. Sousa and P. P. Avelino, Phys. Lett. B689, 145 (2010), 
0911.3902. 

[10] P. P. Avelino, C. J. A. P. Martins, J. Menezes, 
R. Menezes, and J. C. R. E. Oliveira, Phys. Rev. D73, 

123520 (2006), hep-ph/0604250. 

[11] R. A. Battye and A. Moss, Phys. Rev. D74, 023528 

(2006), hep-th/0605057. 
[12] P. P. Avelino, C. J. A. P. Martins, J. Menezes, 



R. Menezes, and J. C. R. E. Oliveira, Phys. Lett. B647, 

63 (2007), astro-ph/0612444. 
[13] P. P. Avelino, C. J. A. P. Martins, J. Menezes, 

R. Menezes, and J. C. R. E. Oliveira, Phys. Rev. D78, 

103508 (2008), 0807.4442. 
[14] P. P. Avelino, J. C. R. E. Oliveira, R. Menezes, and 

J. Menezes, Phys. Lett. B681, 282 (2009), 0907.4389. 
[15] H. Kubotani, Prog. Theor. Phys. 87, 387 (1992). 
[16] B. Carter, Int. J. Theor. Phys. 44, 1729 (2005), hep- 

ph/0412397. 

[17] D. Bazeia, F. A. Brito, and L. Losano, Europhys. Lett. 

D76, 374 (2006), hep-th/0512331. 
[18] P. P. Avelino, R. Menezes, and L. Sousa, Phys. Rev. D79, 

043519 (2009), 0811.1580. 
[19] C. J. A. P. Martins and E. P. S. Shellard, Phys. Rev. 

D65, 043514 (2002). 
[20] T. Boehm and R. Brandenberger, JCAP 0306, 008 

(2003), hep-th/0208188. 
[21] P. P. Avelino, C. J. A. P. Martins, C. Santos, and 

E. P. S. Shellard, Phys. Rev. D68, 123502 (2003), astro- 

ph/0206287. 

[22] P. P. Avelino, C. J. A. P. Martins, C. Santos, and E. P. S. 
Shellard, Phys. Rev. Lett. 89, 271301 (2002), astro- 
ph/0211066. 



